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Abstract 



The elementary excitation spectrum of the spin-i antiferromagnetic (AFM) Heisenberg chain is described in 
terms of a pair of freely propagating spinous. In the case of the Ising-like Heisenberg Hamiltonian spinous can be 
interpreted as domain walls (DWs) separating degenerate ground states. In dimension d > 1, the issue of spinous 
as elementary excitations is still unsettled. In this paper, we study two spin-^ AFM ladder models in which 
the individual chains are described by the Ising-like Heisenberg Hamiltonian. The rung exchange interactions 
are assumed to be pure Ising-type in one case and Ising-like Heisenberg in the other. Using the low-energy 
effective Hamiltonian approach in a perturbative formulation, we show that the spinous are coupled in bound 
pairs. In the first model, the bound pairs are delocalized due to a four-spin ring exchange term in the effective 
Hamiltonian. The appropriate dynamic structure factor is calculated and the associated lineshape is found to be 
almost symmetric in contrast to the Id case. In the case of the second model, the bound pair of spinous lowers 
its kinetic energy by propagating between chains. The results obtained are consistent with recent theoretical 
studies and experimental observations on ladder-like materials. 



1 Introduction 

Low-dimensional quantum antiferromagnets exhibit a rich variety of phenomena indicative of novel ground and 
excited state properties [TJ[51[3]. In one dimension (Id), the ground and low-lying excited states of the spin chain, 
in which nearest neighbour spins of magnitude ^ interact via the antiferromagnetic (AFM) Heisenberg exchange 
interaction, can be determined exactly using the I3ethe Ansatz [H [5] . The ground state has no long range order and 
the spin-spin correlations are characterized by a power-law decay. The elementary excitation is not the conventional 
spin-1 magnon but a pair of spin-^ excitations termed spinous. The physical origin of spinons can be best understood 
in the Ising limit of the exchange interaction Hamiltonian given by 

N ^ J 

H = ^2 J^^i^i+i + 5Z ^ {^tsr+i + s~s^+i) (1) 

i=l 1=1 

where 5,^^ are the spin raising and lowering operators and is the total number of spins. In the Ising limit {Jxy = 0), the 
ground states of the Hamiltonian are the doubly degenerate Neel states, one of which is shown in figure 1(a). An excited state 
is created by flipping a spin from its ground state arrangement, e.g., a down spin is flipped into an up spin in figure 1(b). 
The flip gives rise to two domain walls (DWs) consisting of parallel spins and shown by dotted lines in the flgure Jji6j. The 
transverse exchange interaction term in (1) interchanges the spins in an antiparallel spin pair and has the effect of making 
the DWs propagate independently (figure 1(c)). Since the original spin flip carries spin one, each of the DWs or spinons has 
spin-i associated with it. Spinons are thus examples of fractional excitations in an interacting spin system. 

Spinons can be detected through inelastic neutron scattering in which neutrons scatter against spins to create spin flips. 
Due to energy and momentum conservation, the energy absorption spectrum for spin flips at different wave vectors can be 
measured. One observes a peak at a well-defined energy if the spin fiip creates a single particle excitation. In the case of a 
pair of spinons, the total energy e and momentum k of the spin excitation are given by e{k) = ei(fci) -I- £2(^2) and k = ki + k2 
where ei, ki {i = 1, 2) denote individual spinon energy and momentum ^IjiBj. The total momentum k of the spin flip can be 
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Figure 1: (a) Neel state, (b) a spin flip in the Neel state creates a pair of domain wafls or spinons the locations of 
which are shown by dotted hnes, (c) the transverse exchange interaction term in the Hamiltonian (equation (1)) 
gives rise to the propagation of independent spinons. 

distributed in a continuum of ways among the spinons giving rise to a continuous absorption spectrum. For a single particle 
excitation, the energy versus momentum relation defines a single branch of excitations whereas for spinons a continuum of 
excitations with well-defined lower and upper boundaries is obtained. The compounds CsCoCh and CsCoBrs are good 
examples of Ising-like Heisenberg antiferromagnets in Id above the Neel temperature and provide evidence of the two-spinon 
continuum in neutron scattering experiments [6l[7l[8]. In the case of the isotropic Heisenberg Hamiltonian (Jz = Jxy in 
equation (1)), the spinon spectrum has been clearly observed in the linear chain compound KCuFs [9J though a physical 
interpretation of spinons is not as straightforward as in the Ising-like case. 

The existence of spinons, with fractional quantum number spin-i, is well established in Id Heisenberg- type antiferromag- 
nets. In higher dimensions, the spin-1 magnons are the elementary excitations in magnetically ordered ground states. There 
are theoretical suggestions that quantum antiferromagnets with spin-liquid (no magnetic long range order and without broken 
symmetry) ground states may support elementary spinon-like excitations with fractional quantum numbers [101 111! [12] . A 
well-known example is that of a resonating- valance-bond (RVB) state, a linear superposition of VB states, in which the spins 
are paired in singlet (VB) configurations. A broken VB gives rise to a pair of free spins which may propagate independently 
to give rise to spinon excitations. If the energetic cost of deconfinement is high, the spinons propagate as a bound pair 
(confinement) so that the elementary excitation has spin-1. 

Despite considerable effort, there are few experimental evidences of spinon-like excitations in d > 1 [131 114| . One strong 
candidate is CS2CUCI4, a spin-i Heisenberg antiferromagnet defined on a spatially anisotropic triangular lattice [13| . The 
dynamical structure factor S{k,Lj) for Cs^CuCh, where k and oj are the momentum and energy transfers in the neutron 
scattering experiment, is dominated by a broad continuum which has been cited as evidence for fractionalized excitations. 
Kohno et al [15[ reanalyzed the neutron scattering data to show that the spinons are not characteristic of some exotic 2d 
state but are descendants of the weakly-coupled excitations of individual chains in the material. The spectrum also has a 
sharp dispersing peak attributed to 'triplon' bound states of the spinons. The bound pair lowers its kinetic energy through 
propagation between chains. The issue of fractional versus integer excitations has been extensively investigated in AFM 
spin-i ladders [161 1171 [TSl 1191 120[ . A two-chain ladder consists of two AFM chains coupled by rung exchange interactions. 
The spinons of individual chains are confined even if the rung exchange interaction strength is infinitesimal. Ladders with 
strong rung exchange couplings suppress spinon excitations at all energy scales. Recently, Lake et al [21j have carried out 
neutron scattering experiments on a weakly-coupled ladder material, CaCu203, and shown that deconfined spinons at high 
energies evolve into S = 1 excitations at lower energies. The spinons are associated with individual chains whereas the 
S — 1 excitations are the triplon excitations, i.e., bound states of spinons. Two approaches are usually adopted in probing 
the nature of excitations in spin ladders: (i) the rung coupling strength dominates and (ii) the rung coupling strength is 
weaker than the intra-chain coupling strengths [21|. In this paper, we consider a different case, not studied earlier, in which 
two S = I Ising-like Heisenberg AFM chains, each of which is described by a Hamiltonian of the type shown in equation 
(1), are coupled by Ising or Ising-like Heisenberg AFM exchange interactions. In section II, we investigate the nature of the 
low-lying excitations with Ising rung exchange interactions. In section HI, the rung exchange interactions are considered to 
be Ising-like Heisenberg AFM in nature. Section IV contains concluding remarks. 
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Figure 2: (a) A pair of DWs in the minimum energy configuration of a ladder with an odd number N of rungs, (b) 
Motion of the DW in the top chain leaves in its wake ferromagnetically aligned rung spins which raises the energy 
of the system. 

2 Ising Rung Exchange Interactions 

We consider an antiferromagnetic two-chain spin ladder with the spins of magnitude ^. The individual chains of the ladder 
are described by the Ising-Heisenberg Hamiltonian (equation(l)). The chains are coupled by rungs with the corresponding 
exchange interactions being of the Ising-type. The ladder Hamiltonian Hl is given by 

2 JV JV J 2 N 

Hl = Jz Sj^gSi+i^a + Jz Sf^lSf^2 H 2~'^I/^^ (^ta^i + l,a ~^ ^i,aSi + l,a) (2) 

a — l i=l i=l a — l i — l 

= Hz + HxY 

where the index a — 1(2) refers to the top (bottom) chain of the ladder, i denotes the site index and TV is the total number 
of rungs. We also assume that the anisotropy constant e = ''j^ is << 1. Hence, the Ising part of the Hamiltonian, Hz, 
can be considered to be the unperturbed Hamiltonian with Hxy, containing the transverse exchange interactions, providing 
the perturbation. Since Hz is AFM in nature, the lowest energy states are the Neel states with n.n. spin pairs antiparallel. 
In the spirit of Villain [6], we first consider a ladder with an odd number of rungs, i. e., N = odd and periodic boundary 
conditions (PBCs). Energies are measured w. r. t. that of a Neel configuration of spins. Since there are 3A'' n.n. spin 
pairs, the energy of a configuration in which all such pairs are antiparallel is iSjvee! = — 3A^ Since A*" is odd, a perfect 
Neel configuration is not possible and the lowest energy states of Hz contain a pair of parallel spin pairs which define the 
DWs or spinous (figure 2 (a)). The DWs form a bound pair to ensure minimal energy loss. Any other arrangement of DWs 
in the individual chains gives rise to higher energy states. The lowest energy states are A-fold degenerate as there are A^ 
possibilities for the location of the bound pair which is an Sz = object. 

We next consider the effect of the perturbing Hamiltonian on the minimum energy states. The transverse exchange 
interaction can give rise to independent DW motion in the chains which, however, costs energy as a propagating DW leaves 
in its wake ferromagnetically aligned rung spins (figure 2(b)). The energy cost increases as the distance between the DWs 
increases resulting in confinement of the DW pair. We investigate the dynamics of the DW pair perturbatively using a low- 
energy effective Hamiltonian (LEH) |22II23| . The A-fold degenerate DW pair states \pi}, i ~ 1, N (figure 2(a)) constitute 
the low-energy manifold and have energy Eg = — ^'^^^ + Jz- The higher energy states of Hz are denoted by \qc,) with energy 
Ea- The perturbing Hamiltonian Hxy connects the low-energy manifold to the manifold of higher-energy states. In general, 
perturbation lifts the degeneracy of the low-energy manifold leading to an effective Hamiltonian operating in the space of 
states associated with the low-energy manifold. Diagonalization of the n-th order (n — 1,2, ...) effective Hamiltonian in the 
low-energy subspace of states reproduces the n-th order energy corrections to the low-energy unperturbed states. Using 
degenerate perturbation theory, the first order LEH is given, up to an overall constant, by [23j 
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Figure 3: Two successive applications of Hxy on low-energy states shifts the location of the bound DW pair by 
two lattice constants. 



The second-order LEH has the form 

ij a 

Since the matrix element {pi\HxY\pj} ~ 0, the LEH is determined using the second-order expression in equation (4). Two 
types of processes contribute to H^^j:. In "diagonal" processes, the spins in an antiparallel pair exchange and then reexchange 
back to the original configuration {\pi) = \pj))- Such processes do not lift the degeneracy and give rise to a constant energy 
shift. We neglect this contribution in deriving the effective Hamiltonian. In the off-diagonal processes, the spins in two 
antiparallel pairs belonging to the four-spin plaquettes bordering the bound DW pair are interchanged. For the DW pair 
state shown in figure 2(a), the intermediate states \qi), i = 1, .., 4 are given by 

ln\ t ^ 4 t i t ; f.. 



k2) 



t i t t i t i 

i 4 ^ i t i t 

t i t 4 ^ t i 

i t i i t i i 

t i t t i t i 

i t i ^ 4 i t 



The spin pairs deviated from the arrangement shown in figure 2 (a) are marked by double arrows. The perturbing Hamilto- 
nian, Hxy, acting on the intermediate \qa) states shifts the location of the bound DW pair by two lattice constants either 
towards the left or the right. This is illustrated in figure 3 for the state \qi)- The energy of the states \qa) {a = 1, ...,4) is 
-Eq = -^^^ + 2Jz. The second-order LEH is thus given by 

Hiff = E i\P^+^)(P^\ + b-2>(P»l) (6) 

i 

where e = ''j^ . The off-diagonal processes are equivalent to "ring" exchanges involving four spins. The full second-order 
Hamiltonian, defined in the low-energy manifold, is thus given by 

Heff — Hz + Hring (7) 

where 

Hring ~ Jring ^ ^ ('^1 '^2 "^3 "^4 h.C.^ (8) 

□ 

with Jring = and the sum over all elementary plaquettes of the ladder. Ring or cyclic exchange interactions 

(equation (8)) also appear in the perturbative effective Hamiltonian theories developed for the XXZ Heisenberg model on 
the checkerboard lattice [24] and in the case of an easy-ajds Kagome antiferromagnet [25| . In the ladder model, the ring 
exchange interaction has the effect of deconfining the bound DW pair. In the low-energy subspace, the dispersion relation of 
the bound pair can be determined in a straightforward manner. The eigenstate tp{k) of the pair can be written as a linear 
combination of staes \ j){j = 1, 2, A'^) where j denotes the location of the bound DW pair. 

1 

m)} = ^J2^"''\j) (9) 
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(c) 

Figure 4: Bound DW pair states for (a) fi — 1 and (b) /i = 3. The dotted lines indicate the locations of the bound 
DW pairs; (c) One of the ground states, tpNeei i, of Hz- 

Heff (equation (7)) operating on \^{k)} yields the eigenvalue 

ujb{k) = J z{l ~ cos 2k) (10) 

The subscript '6' in uj^ denotes that the dispersion relation is that of a bound DW pair. 

We next consider a two-chain spin ladder with an even number A'' of rungs and described by the Hamiltonian in equation 
(7) satisfying PBCs. The low-lying excitation spectrum is obtained in the subspace of degenerate eigenstates of Hz which 
are generated by flipping all the spins in a block of ^ (/i may be odd/even) adjacent rungs in the ground state (Neel state) 
of Hz- Each of the states contains two bound DW pairs (figure 4) and has energy 

Sdvk = -^^ + 2Jz (11) 
4 

The z-component of the total spin of each state, 5*°' = 0. We consider /i to be odd with the degenerate eigenstates of Hz 
given by 



V'i(fc) = 

Mk) = 
^pN-i{k) = 




where Sk ~ S^iSj. 2 and 5^ — S/. "fpNeeii is the ground state of Hz shown in figure 4(c). The Hamiltonian Hring in 

(7) has the following matrix elements between the Ising eigenstates : 

{Mk)\H„r.a\Mk)) = {Mk)\H„r.3\^pi{k)) = ... = J™^ (l + e""'") = v (13) 
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Figure 5: Excitation continuum with energies given in equation (17) and e = 0.15 



where Jring 



2 • 



The low-lying excited state of the Hamiltonian H^ff (equation (7)) is given by 

N/2 

1pDw{k) = Ci.^2i/-l(fc) 



(14) 



From the eigenvalue equation HefftpDw{k) = \^Dw{k), one obtains 

N/2 



(15) 



where 

{u\H,ff\u') = 2Jz foru' = v (16) 
= V for i/' = u + 1 
= V* for v' = V — 1 
= otherwise 

The diagonal matrix element is 2Jz as energies are measured with respect to the Neel state energy — '^^^^ (see equation 

The eigenvalues constitute an excitation continuum given by 

■ cos A: cos(fc + 0)1 (17) 



(11)). We choose the coefficients C^s to be Cv = e 

\ = 2Jz 



where — tt < (p < n. 

Figure 5 shows the excitation continuum with upper and lower bounds given by 27^(1 ± e^cosfc). The degenerate 
eigenstates of Hz defined in (12) correspond to the top chain of the spin ladder being in as Stot,i = +1 and the bottom 
chain being in an Stot,2 = — 1 state. One can construct a set of degenerate eigenstates with the situation reversed. Also, /x, 
the number of adjacent rungs constituting the block of fiipped spins can be even {fj, = 2, 4, 6, ...etc.). There are two distinct 
sets of such states |8l[2^. All these subspaces of states give rise to the same excitation continuum (figure 5). The excitation 
continuum arises due to the motion of two walls each of which consists of a bound pair of DWs. 

One notes that the effects of Bring, for the two-chain ladder, and Hxy, for the Id chain, on the DW states are similar. 
In the first case, a bound pair of DWs shifts by two lattice constants and in the second case a single DW shifts by the same 
distance. In the case of a single chain, the eigenvalues Aid constituting the excitation continuum are 



Aid = Jz [1 + 2e cos k cos(fc 4- ( 



(18) 



Comparing equations (17) and (18), one finds that the spread of the continuum around the unperturbed level is less in the 
case of the spin ladder. 

The dynamic form factor, Si\{q,u]), associated with the bound DW pair is defined at T = by 



Sii(g,a;) = ^ \{f\A{q)\g)\' 5{i^ - Ef + Eg 



(19) 



where 



.(.)^^E-^^---(^^^^ + ^a^^) 



(20) 
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In (19), If?) and |/) are the ground and excited states of H^ff connected by A{q), with energies Eg and Ef respectively, uj 
and q are the frequency and wave number of the excitation. Sii{q,ijj), involving a pair of spin deviations, could be probed 
by the hght-scattering techniques [ST] [28]. Upto the first order of Jring [Jring << Jz in (7)), 

\g) ~ tpNdel 1 + T^HringlpNeell (21) 

EjQ — tlz 

where Eo is the energy of ipNeeii- Since Hring acting on ipNieii creates two bound DW pairs, g^j';^^^ ~ — jTJ-Thus, 

A{qM ~ -J= (^1 + ^ COSq^ Ml) + -^^^^^iq) (22) 

where tpi{q) and tp3{q) are as defined in equation (12) and v* = — ^ (l + e^^'^) . Using equation (22) and (19) and the 
expression (14) for |/) — tpDw{k), one gets following the procedures described in [71 126] 

S^,{q,u;) ^ ^l^Z^(^l + e'cosq--^^for\Q\<2\v\ 

= otherwise (23) 

with = IX) — 2Jz ■ 

The expression (23) is similar to that for the dynamic structure factor S^^iq,^) of the Id chain obtained in first order 
perturbation theory [7] {Ai = in equation (20)) except that in the latter case, n — uj ~ Jz and the contribution of the 
anisotropy term is to first order in e. Figure 6 shows the plots of Sii{q,uj) x 2Jz \ cosq\ versus -j^ for e = 0.15 and for various 
values of the wave number q. The lineshape is almost symmetric in contrast to the prominent asymmetry found in the Id 
case [7]. 



3 Ising-Heisenberg Rung Exchange Interactions 

We now consider the case in which the rung exchange interactions of the two-chain spin ladder are Ising-like Heisenberg-type. 
The Hamiltonian is given by 

2 N N J 2 N 

Hr = Jz'^'^SlaSi + l,a + Jz'^SliSl2 ^ ^"^"^ {St^aSi + l^a + Si^aSt+l,c,) (24) 

a—l i—l i—l a—l i—l 

JXY ^ 

i = l 

= Hz + HxY 

The Hamiltonian (24) differs from Hl is equation (2) by the addition of the last term. The ground states of Hz are the 
doubly degenerate Neel states. We consider the Neel state iliNieii shown in figure 4(c). The ladder can be divided into 
two sublattices A and B such that in tpNeeti the A{B) sublattice spins are pointing up (down). The ground state energy 



3NJz 
4 

The lowest energy excitation of the unperturbed Hamiltonian is obtained by flipping a single spin in either the A {Stot 



Eo = 4 



— 1) or the B {Stot = +1) sublattice. We consider the latter case with \i) denoting the state in which the flipped spin is 
located in the ith rung (figure 7(a)). These excited states are A'^-fold degenerate with the energy 

+ ^ (25) 

The perturbing Hamiltonian acting on the state |i) generates the following states 

HxY\i}='-^[\l} + \2} + \3} + \4} + ....] (26) 

where e = '^fj . The states |m)(m = 1, 4) are shown in figure 7(b) with energy 

_ 3NJz 5Jz , „s 

iim h — (27) 

The other states which are generated when Hxy acts on the state \i) have higher energies and are hence not considered. 
HxY acting on the states |m) gives 

HxY\m) = l:^(li)+lj-l)), m = l,2 (28) 
= ^(l^) + |^ + l)), m = 3,4 
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Figure 6: Sii{q,uj) x 2Jz\ cosg| evaluated in second-order perturbation theory for e = 0.15 versus for different 
values of q. 
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In first order perturbation theory there is no energy correction. A finite energy correction is obtained in the second order 
perturbation tlieory. The effective LEH (equation (4)) with \pi} = \i}, \pj} = \j}, \qa} = |m) and Ea = Em in the same 
order is given by 



2 

2 

i 

N 



H% = + + (29) 



~ 2 — 51/ 51/ iStl^i+S,2 + Si,lSt+S,2) (30) 
i=l (5=-l,l 

The effect of this Hamiltonian on the low-energy excited state |i) (figure 7(a)) is to shift the flipped spin from the ith to the 
(i + l)th or the {i — l)th rungs. Since the flipped spins are located in the B sublattice, the shift is in the diagonal direction. 
The flipped spin is associated with a bound pair of DWs in a chain. The bound pair lowers its kinetic energy by propagating 
between chains. The full second-order Hamiltonian defined in the low-energy manifold of states with single spin fiips, is 

H.ffi =Hz + H^ff, (31) 

As before, we have not included the terms arising from the "diagonal" processes in equation (4) as they give rise to a constant 
energy shift. The low energy excited state with Sfot = +1 can be constructed as 

1 ^ 

where 'B' denotes the B sublattice. The dispersion relation for the propagation of the flipped spin or equivalently the bound 
DW pair is given by 

Ei{k) = Jz {I- cos k) (33) 

where the energy is measured w. r. t the Neel state energy. The bound DW pair moves diagonally across the spin ladder. 
The unperturbed Hamiltonian, Hz, is the same irrespective of whether the rung exchange interactions are Ising-type or 
Ising-Heisenberg-type. Thus, the lowest unperturbed excited state is the single flip state in both the cases. The excitation 
has a localized character when the rung exchange interactions are of the Ising-type. The bound DW pair associated with 
the single spin-flip can not propagate between the chains as the inter-chain interactions are Ising-like and propagation of the 
DWs in a single chain is, as pointed out before, energetically prohibitive. The energy, — of the localized excitation is 
lower than that of the propagating excitations involving two bound DW pairs when the rung exchange interactions are of 
the Ising-type (section 2). In this case, propagating excitations with the lowest energy involve two bound DW pairs rather 
than one. 



4 Concluding Remarks 

AFM spin models in which the existence of spinous is well-established include the spin-i Heisenberg AFM chain [1], the 
Majumdar- Ghosh model [29] and the Haldane-Shastry model [301 131| . The physical picture of a spinon as a DW between 
two degenerate ground states emerges in the Ising-Heisenberg limit of the AFM Hamiltonian |T][7l[8|. In the case of the MG 
model, the spin-^ excitation acts as a DW between the two dimerized ground states of the model [321 133] . In a closed chain, 
the DWs occur in pairs so that the lowest-lying excitation is given by the two-spinon continuum. The spinous are deconfined 
in this case and can move away from each other. There is no energy cost in moving the spinous far apart. This is not so when 
two AFM chains are coupled in the form of a spin ladder. Let Jxand J\\ be the strengths of the rung and intra-chain n.n. 
exchange interactions respectively. The spinon excitations of individual chains are confined by even an infinitesimal coupling 
strength Ji |18ll2T] . The two 5 = ^ spinous form a bound state giving rise to singlet and triplet excitation branches. In the 
case of the strongly coupled ladder {J± >> J||), the elementary excitation is a triplet. Lake et al. |21) carried out neutron 
scattering experiments on the weakly-coupled (J_l << J\\) ladder material CaCu203 and obtained evidence of the singlet 
excitation mode. The spinon continuum was observed at high energies for which the chains are effectively decoupled. The 
spinous in a chain evolve into an S = 1 excitation at lower energies thereby confirming that the S = 1 "triplon" excitation is 
a bound state of two spinous and not a conventional magnon. 

In this paper, we study a two-chain S — ^ AFM spin ladder in which the individual chains are described by the Ising-like 
Heisenberg Hamiltonian and the rung couplings are of the Ising-type. Using a low-energy effective Hamiltonian approach, 
we establish that in a ladder with an odd number of rungs the spinous (DWs) form a bound pair. A four-spin ring exchange 
interaction in the effective Hamiltonian is responsible for the delocalization of the bound pair. In the case of a ladder with 
an even number of rungs, the low-lying propagating excitation involves two bound pairs of DWs which can move away from 
each other giving rise to a continuum of excitations. The physical origin of the excitation continuum is similar to that in the 
case of the Ising-Heisenberg AFM chain in Id except that in the former case the spinous form bound pairs and the dispersion 
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Figure 7: (a) The lowest energy excitation {S^^^ = +1) of the unperturbed Hamiltonian Hz in equation (24). 
The deviated spin is represented by a thick arrow, (b) The perturbing Hamiltonian, Hxy, acting on the state \i) 
generates the states |m)(m — 1,2,3,4). The spin deviations from the state |i) are shown by thick arrows. 
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of the excitation spectrum is a higher-order effect in perturbation theory. This results in an almost symmetric lineshape in 
the case of the dynamic structure factor Sii{q,u}) (Figure 6) in contrast to the asymmetry observed in the structure factor 
Sxx{q,oj) in Id 0- The delocalization of a bound spinon pair is brought about via a ring or a diagonal exchange interaction 
term in the effective Hamiltonian. 

We further consider a second model in which the rung exchange interactions are described by the Ising-like Heisenberg 
Hamiltonian. In this case also, the spinon pair in a single chain is bound and the bound pair lowers its kinetic energy by 
hopping between chains. Kohno et al. [15| have studied a, S — ^ spatially anisotropic frustrated Heisenberg antiferromagnet 
in 2d in the weak interchain coupling regime. The model provides a good quantitative fit to the inelastic inelastic neutron 
scattering data of the triangular antiferromagnet Cs2CuCh. The spectrum consists of a continuum arising from the decon- 
finement of spinous in individual chains and a sharp dispersing peak associated with the coherent propagation of a triplon 
bound state of two spinons between neighbouring chains. In the case of our model, the bound pair has = +1. One 
can similarly construct an = — 1 excitation. In summary, we have studied ladder models with the Ising-like Heisenberg 
Hamiltonian describing the interactions in individual chains. The rung interactions may be pure Ising or Ising-like Heisen- 
berg. The ladder models studied in this paper share some common features with models in which the exchange interactions 
are isotropic. In the latter case, two types of models have generally been considered : ladder models in which the rung 
exchange interactions are the most dominant and models which describe spin chains coupled by weak exchange interactions. 
The models considered in this paper belong to a category not studied earlier and provide considerable physical insight on 
the origin of spinon confinement and how the bound spinon pairs delocalize. The isotropic models have, however, a richer 
dynamics with interactions generating cascades of virtual particles so that the two-body confinement problem becomes a 
many-body one |21) . The major common feature emerging from the study of ladder models with both Ising-like and isotropic 
exchange interactions appears to be the confinement of spinons in the form of bound states. The origin of excitation continua 
in specific cases lies in multi-triplon excitations rather than in the fractionalization of excitations [18| . 
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